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1. Motivation and objectives 

1.1 Introduction 

Although the different regimes of premixed combustion are not well defined, most 
of the recent developments in turbulent combustion modeling are led in the so- 
called flamelet regime. The goal of these models is to give a realistic expression to 
the mean reaction rate ( w ). Several methods can be used to estimate (xb). Bray 
and coworkers (Libby &: Bray 1980, Bray 1985, Bray & Libby 1986) express the 
instantaneous reaction rate by means of a flamelet library and a frequency which 
describes the local interaction between the laminar flamelets and the turbulent flow 
field. In another way, the mean reaction rate can be directly connected to the flame 
surface density (E). This quantity can be given by the transport equation of the 
coherent flame model initially proposed by Marble & Broadwell 1977 and developed 
elsewhere (Pope 1988, Candel Sz Poinsot 1990, Trouve & Poinsot 1992, Bidaux et ai 
1993). The mean reaction rate, (u;), can also be estimated thanks to the evolution 
of an arbitrary scalar field G(x, t) = Go which represents the flame sheet. G(x, t) is 
obtained from the G— equation proposed by Williams 1985, Kerstein et al. 1988 and 
Peters 1993. Another possibility proposed in a recent study by Mantel & Borghi 
1991, where a transport equation for the mean dissipation rate ( e c ) of the progress 
variable c is used to determine (w). 

In their model, Mantel &; Borghi 1991 considered a medium with constant den- 
sity and constant diffusivity in the determination of the transport equation for (c c ). 
A comparison of different flamelet models made by Duclos et al 1993 shows the 
realistic behavior of this model even in the case of constant density. Our objective 
in this present report is to present preliminary results on the study of this equation 
in the case of variable density and variable diffusivity. Assumptions of constant 
pressure and a Lewis number equal to unity allow us to significantly simplify the 
equation. A systematic order of magnitude analysis based on adequate scale rela- 
tions is performed on each term of the equation. As in the case of constant density 
and constant diffusivity, the effects of stretching of the scalar field by the turbu- 
lent strain field, of local curvature, and of chemical reactions are predominant. In 
this preliminary work, we suggest closure models for certain terms, which will be 
validated after comparisons with DNS data. 

1.2 The role of scalar dissipation for premixed turbulent combustion 

For the prediction of chemical reactions in premixed media, the dissipation time 
scale of the concentration (or the temperature) fluctuations r c and of the turbulent 
kinetic energy r% are crucial quantities. Nevertheless, in most turbulent combustion 
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models (flamelets models, presumed pdf approach) where these scales are taken into 
account, proportionality between these quantities is very often assumed. Although, 
numerous experimental studies show that the ratio of these two scales is not a 
constant and varies for different flows (see Beguier et al. 1978). 

One way of calculating the dissipation time scale of a scalar is to use a transport 

dc* dc * 

equation for the dissipation of the fluctuations of this scalar (e c ) = 

(where d is the molecular diffusivity and c' the fluctuation of c; c is defined to be 
equal to 1 in the burnt gases and to 0 in the fresh gases). Starting from the studies 
of Zeman & Lumley 1976, Borghi & Dutoya 1978 used this equation to investigate 
the case of a steady premixed flame. But the use of such am equation requires the 
modeling of many terms and very little knowledge exists to assess the adequacy of 
these models. 

More recently, Mantel & Borghi 1991 have developed a transport equation for 
(e c ) in the case of constant density and constant diffusivity. An order of magnitude 
analysis based on appropriate scale relations shows that the effects of stretching 
of the concentration field by the small scales of turbulence and the effects of local 
curvature are predominant. The authors present new closures for these terms in 
order to take into account the presence of laminar flamelet within the flow. 

The applications of such an equation are significant and of fundamental interest: 

1- In the particular case of high Damkohler number, Da = T|/r c fc where r c /, is a 
chemical time scale, Bray & Moss 1977 and Borghi 1978 have demonstrated that 
( w ) and (e c ) are directly related: 


(w) = 


(ec) 

[~bo 


( 1 . 1 ) 


where, 


bo — 



In accordance with Bray Si Moss 1977, the equation for (e c ) leads to a transport 
equation for the mean reaction rate and to a new formulation of the well known 
Eddy Break-Up model. 

2- Moreover, in the flamelet regime, (e c ) and (E) can be directly connected by 
(e c ) oc 5^(E) (see Borghi 1990), where Sl is the laminar flame velocity. Thus, 
a transport equation for (E) is obtained in a different way from those proposed 
by other studies (Marble Si Broad well 1977, Candel Si Poinsot 1990, Trouve Si 
Poinsot 1992, Bidaux et al. 1993). 

3- Furthermore, the knowledge of the dissipation time scale for the scalar allows 

Mantel et al. 1993 to propose a new model for the equation for the correla- 
tion which is another crucial problem in turbulent modeling that is often 

underestimated . 
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2. Accomplishments 

We present m this section the different steps and assumptions leading to the 
equation for (e c ). 

The methodology chosen here to calculate the transport equation for (E c ) is 

to first establish the instantaneous equation for E c = d— — . An averaging 

operation and an order of magnitude analysis of this equation will lead directly to 
the equation for ( e c ). 


2.1 Preliminaries 

In order to determine the transport equation for E c , some assumptions are made. 
A perfectly premixed medium is considered and a single chemical reaction occurs 
m the flow which can be described by: F uel + 0 xydlzer - P roduct8 . 

From the instantaneous transport equation for the progress variable c: 


— + u dc - 1 

dt * dxk pdxk 




dx k 


( 2 . 1 ) 


The transport equation for E c is obtained by differentiating Eq. (2.1) with respect 
to x, and by multiplying the result by dc, „ where c, represents After some 

algebra, the transport equation for E c leads to: ' 


dE c 

dt 


+ U k 


dEc 

dx k 


= + U k d tk c <iC<i - 2dc t ,c }k ~ + -^-( p d^k) 

oxi poxk' dx */ 

~ d k d k c,c t , + 2dd, ,c~ -2 d 2 p-p-+ 2 dc ,w , 

ax k ax k ' ’ 


,dE c 


'dx\ 


Obviously, Eq. (2.2) involves unclear correlations and cannot be used directly. New 
assumptions are needed in order to simplify Eq. (2.2) and to clarify the physical 
meaning of its terms. Let us consider setting the Lewis number equal to unity, 
no heat losses, and a constant pressure. Concerning this last assumption, even if a 
pressure drop occurs across a laminar flame, it is so weak that the density variations 
are essentially due to temperature variation across the flame front. This assumption 
could only be put into question in case of strong pressure waves spreading into the 
fiame front (i.e. knock in I.C. engines). 

In case of a Lewis number equal to unity and no heat losses, the temperature and 
the progress variable are directly related by: 


T = T o (l + 0c) (2.3) 

where 8 is the heat release parameter: 8 = T ad and T 0 represent respec- 

tively the adiabatic temperature and the temperature of the fresh gases. 
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Due to assumptions of constant pressure, unit Lewis number, and no heat losses, 
density and diffusivity derivatives can be related to the progress variable derivatives 
which allow to simplify Eq. (2.2). 

The assumption of constant pressure and Eq. (2.3) imply: 



2 P 


P,k = ~P 

— 

(1 +0C)2 


9 

P 1 + 9c C ’ k 

(2.4) 

e dc t i 

2 c,c* P l+ 0 cdx k 

(2.5) 


and allow us to replace in Eq. (2.2) the density derivatives by the progress variable 
derivatives. 

The diffusivity derivatives can be solved by considering that the dynamic viscosity 
H is only a function of the temperature, for example: 


P - Po 



( 2 . 6 ) 


where b is a constant equal to 0.76. The dynamic viscosity at temperature T 0 is 
denoted by po- Thus, the first and second derivatives for the molecular diffusivity 
can also be related to the first and second derivatives of the progress variable: 


d ‘ = ll + k}d TTii c ’ 1 

S = ^ + 1) (TW dw + ( ‘ + 1) T^£i 




0 2 0 dc i 

(TTW dc ' iC ' k + { b + l) TT^cdr k 


(2.7) 

( 2 . 8 ) 
(2.9) 


Substituting these new expressions for the derivatives of the density and molecular 
diffusivity, Eq. (2.2) becomes: 


5E c , . 3E C d ( , dE c \ j dUk ^ 

w + pUt W k “ dT t y pi 9T t ) ■ ipic ‘ c ' k + (1 + )p ‘ 


dx k 

2bp 


dUk 

dxk 


9 2 


(1 + 9c? 


„ 9 d 2 c 

Ec + {l + b)p TT9~c dEc fa 2 . 


2(2 + b)p 


9 ,o dc i 2 ^ C .» ^ c .« , O j 

d c,i c ,kirr - 2 P<p7rr*± + 2 P dc ,> w P 


1 0c 


dx k 


dx k dxk 


( 2 . 10 ) 


In Eq. (2.10), derivatives of density and diffusivity disappear and a physical 
meaning of the terms of Eq. (2.10) can be seen more clearly. But before we give 
physical meanings to these terms and propose closures, averaging operations and 
an order of magnitude analysis are performed. We will see that the instantaneous 
gradients of c are predominant and allow us to derive an equation for (e c ). 
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2.2 Averaging and order of magnitude analysis 
2.2.1 Preliminaries and scale relations 

First, we decompose the velocity and the progress variable into a mean value 
and a fluctuating component (i.e. for c, we have c = (c) + c ' where (c) = ^ is 

the Favre averaging, and by definition (c') = 0). The order of magnitude analysis 
needs appropriate scale relations which take into account the presence of chemical 
reactions in the fluid. For the gradient of the velocity, the approach of Tennekes Sz 
Lumley 1972 (pp 88-92) for the vorticity equation is retained. These scale relations 

A/rr \ t a # . 


d(Uk) u' j du k u' 

3X6 q x . 0C Y t 3X1(1 x wllere ‘t and A are respectively the integral length 

scale and the Taylor microscale of the turbulence. 

For the quantities related to the concentration flow field, the problem is more 
subtle. Because chemical reactions occur within the flow, we have to take into 
account the existence of flamelets in the turbulent flow field in order to estimate 
the scale relations, and we must also know the pdf (probability density function) of 
c for the averaging operation. 

In the limit of very large Da, this pdf is essentially composed of two peaks cor- 
responding to the fresh and burnt gases: 


P{c) = S 0 S(c) + 7 /(c) + <M(1 - c) 


( 2 . 11 ) 


The constant 7 in Eq. (2.11) must to be proportional to Da -1 as demonstrated by 
Bray & Moss 1977. This can be easily justified by noting that in the case of very 
fast chemical reactions (r c — > 0), the instantaneous reaction rate tends to infinity 

{ T 7h = / w(c)dc) while the mean reaction rate, defined by: 

J 0 

(w) = f w(c)P(c)dc (2.12) 

Jo 

has to remain finite (since u>(0) = 0, w( 1) = 0). This point has a strong impact on 
the order of magnitude analysis on terms where concentration gradients appear. 
We now go back to the gradients of the progress variable. Mantel and Borghi 

1991 proposed to relate mean gradients to the integral length scale oc 1 

dxi l t 

fk f / t2 \ 1 /2 

and instantaneous gradients to the laminar flame thickness 61: — oc — — - 

dxi 61 

The relations Si oc (d/r^) 1 / 2 and Si oc ( dr c * ) 1 ! 2 are also needed in this analysis. 
As an example, we can examine the case of (E c ) itself. By definition, we have: 


(Ec) = [ dc,,c t P(c)dc = d^Q-Da- 1 = ^l D a 0 Re° t 


n 


(2.13) 


where, Re t is the Reynolds number based on the integral length scale 
Re t = u'lt/u. 
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Decomposing c into a mean value and a fluctuating component, (E c ) can be 
developed in different terms: 

(E c ) = (dc t ic t i) — (d)(c > i)(c i j) + 2 (dC'j)(c t i) + (dc ti c ti ) ^ 

Re~ l Da 0 Re~ l, 3 Da~ l/i Re®Da® 

and we can conclude that in the case of large Reynolds and Damkohler numbers, 
only the last term of Eq. (2.14) has to be retained: 

(E c ) = (dc'iC,) = (e c ) (2.15) 

This explains why we can derive an equation for (e c ) from the equation for (E c ). 
2.2.2 The averaged equation 

The separation of the velocity Uk and the progress variable c into a mean value 
and a fluctuation leads to a large number of correlations and new terms in the 
Eq. (2.10). However, the order of magnitude analysis permits an investigation 
concerning the relative importance of the terms of the equation and allows removal 
of these terms which are negligible. 

After a systematic study of each term in Eq. (2.10), we obtain: 




dt 


dx k 


(i) 


(II) Re°Da° 


dx k 

(III) Re~ U 1 Da 1/2 




(IV) Re® Do 0 


+ bp{€c ^ ~ P ^ dC ' iC ’ k)J ^)* k P ( {dC ’ iC ' k dx ‘0 

(VIII) Re, I/l Da' 

dc,i€ c ) 


dx k 

(V) Re® Da® 
i^k 


(VI) Re} / 3 Da ( 
a 2 

- 2 bp( 


(l + Sc)* 
(IX) Re®Da 


(VII) Re“Da° 
e 2 c ) - 4bp(c,)( 




(1 + Bcf 

(X) Re," 1 /J Da l/J 


a pa r > Q Q 2 C » 

- 2(6 + 3)?(c,}( rr ^,iV, g - f > - (6 + 3 

(XI) Rt; 1 , 2 Da 1 ' 7 (XII) Re® Da 1 

,-,ji 3c,i dc t i 


- 2 p(d 2 


dxk dxk 

(XIII) Re® Da 1 


) +2 p(C'i){dw ti ) +2p{c[ i dw[ i ) 

(XIV) Rt~' n Da l/J (XV) Re 0 , Da 1 


(2.16) 


Some physical interpretations are proposed for the terms appearing in Eq. (2.16). 
From left to right, the following terms appear: 

- (I) accumulation 

- (II) convection 
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- (Ill) molecular diffusion 

- (IV) turbulent diffusion 

- (V) production (or destruction) due to expansion (or compression) of the mean 
velocity field 

- (VI) production (or destruction) due to expansion (or compression) of the turbu- 
lent flow field 

- (VII) production due to mean velocity gradient 

- (VIII) production by stretching of the concentration field by the turbulent strain 
field 

- (IX) destruction 

- (X) production (or destruction) due to mean concentration gradients 

- (XI) production (or destruction) due to local curvature effects 

- (XII) destruction due to local curvature effects 

- (XIII) destruction due to molecular dissipation 

- (XIV) production (or destruction) due to chemical reactions 

- (XV) production (or destruction) due to chemical reactions 

2.8 First closure attempts 

In this section, we present closures for some of the terms of Eq. (2.16). Because 
this study is in a preliminary phase, closure models do not exist for some terms but 
are currently under investigation. 

- The turbulent diffusion term (IV) is closed by relating the correlation {u' k e c ) to 
a turbulent diffusion coefficient d t and to the gradient of (e c ): 


— ( u 'k e c) 


dt d(t c ) 

<7 C dx k 


(2.17) 


where a c is a Schmidt number. 

- The term due to expansion or compression (VI) can not be closed easily. The 
sign of this correlation depends on the local velocity field. This sign could be 
determined by analysis of DNS results. 

- The production term due to mean velocity gradient (VII) can be modeled as in 
Zeman & Lumley 1976: 


~P 




(“X) d(Uk) 

k dxi 


(2.18) 


where Cp v =1.0 and k represents the turbulent kinetic energy. 

- The production by stretching of the concentration field due to the turbulent strain 
field (VIII) can be modeled as in the case with constant density. Because the sign 
of this correlation seems to be negative (see the experimental results of Antonia 
& Browne 1983), this term represents a production of ( e c ). Thus, the model 
proposed by Mantel & Borghi 1991 is still valid: 


-P 



i^k 


= a 0 pRe] /2 ^- 

Tt 


(2.19) 
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where Qo is a constant (oq = 0.9). 

Q2 

- The term —2bp(y- 1S al wa y s negative and represents a destruction of 

( 1 + vC) 

(e c ). We propose as a first approximation: 


2 bp( 


0 2 


(1 +0cf 


t\) = R lP 


e 2 


(l + 9(c)) 2 


< £ c ) 2 


( 2 . 20 ) 


where R\ is a constant. Its value could be estimated from DNS results. 

- The production (or destruction) due to mean concentration gradients 

02 

4 bp(c t i)(— — — dc^e c ) can not be easily modeled. More information concerning 

(1 + vcy 

the correlation c % {€ c has to be extracted from DNS or analytical study. 

_ 0 <9V 

- The term —2(6 + 3)p(c i ,)( ^ * s a P roc ^ uct * on ( or a destruction) of 

(e c ). No clear closure can be proposed, and more investigation must be pursued. 

- The term (b + 3)p( - de c is a production (or destruction) of (e c ) due to 

local curvature effects. Following the closure proposed by Mantel & Borghi 1991 
for the curvature, we propose: 


0 


& 2 


0 


( £ c ) 2 


(& + %( 1 + 0 c * c 5x 2 > R2p i + e( c ) {c>*y/* Retl ^ 2 - 21 ^ 


d 2 c' 

R 2 is a constant. Its sign depends of the sign of the correlation e c -^— 2 ” itself which 

ox ^ 

could be obtained using DNS results. 

For the term of destruction due to molecular dissipation (XIII), the model pro- 
posed by Mantel &: Borghi 1991 is still valid in the case of variable density. From 
these authors, we have: 


2 p{d 2 


dc )t dc )X 
dxk dxk 


= fiopRe 


1/2 (ec ) 2 

1 (C 2 ) 


( 2 . 22 ) 


fio is a constant (/?q = 1.25). 

- concerning the production (or destruction) terms due to chemical reactions (XIV) 
and (XV), further investigations are in progress in order to propose new closures 
for these terms. 


3. Conclusion and future plans 

A transport equation for the dissipation of the fluctuation of a reactive scalar 
has been established in the case of variable density and variable diffusivity. The 
assumptions of constant pressure, no heat losses, and a Lewis number of unity 
significantly simplify the equation. As in the case with constant density (Mantel &; 
Borghi 1991), an order of magnitude analysis using adequate scale relations shows 
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that the effects of stretching, molecular dissipation, local curvature, and chemical 
reactions are dominant. As a first approximation, some closure models are proposed, 
but further theoretical developments and comparison with DNS results must be done 
to confirm these models. 

Analytical developments must be in the limit of large Damkohler and Reynolds 
numbers in order to study the relative importance of the unclosed terms. 

DNS could also give essential information on the sign of some correlations and 
will allow the study the different terms of the equation separately. Existing DNS 
results of flame propagation in a decaying homogeneous turbulence for different 
Lewis number will be used. New computations will be performed over a large range 
of the ratio if' /Sl in order to assess the validity of the models and to complete the 
CTR DNS database. 
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